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Abstract 

Riemannian metric on real 2n-dimensional space associated with 
the equation governing complex diffusion of pure states of an open 
quantum system is introduced and studied. Examples of a qubit under 
the influence of dephasing and thermal environments are used to show 
that the curvature of the diffusion metric is a good indicator of the 
properties of the environment dominated evolution and its stability. 
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1 Introduction 



The states of an open quantum system are commonly described by a density 
matrix p. In many cases, the evolution of p(t) is governed by a master 
equation of the Linblad form p],[2] for the density matrix p(t) 

^ = -i[H,p] + \Y\Uf>iL\\ + [Li,0ii (1) 

where the Linblad operators L\ describe the influence of the environment. 
The equation (1) represents the general form of an evolution equation for a 
quantum system which satisfies Markov property. 

However, this theoretical approach to the dynamics of open quantum 
systems is not unique. In real experiments it is often useful to understand and 
model the dynamics of pure quantum states [3], [I], [5]. Indeed, the evolution 
of an open system can be described directly in terms of the dynamics of the 
system's pure state. The corresponding evolution equation is a stochastic 
modification of the unitary Schroedinger equation. In fact, the density matrix 
p can be written, in different but equivalent ways, as a convex combination 
of pure states. Each of these results in a stochastic differential equation for 
\ip{t) > in the Hilbert space TC. Such stochastic Schroedinger equations (SSE) 
are called stochastic unravelling [6], [7], [2] of the Linblad master equation for 
the reduced density matrix p(t). There are many different forms of nonlinear 
and linear SSE that have been used in the context of open systems [E],[3], 
[2j , [6] , [7] or suggested as fundamental modifications of the Schroedinger 
equation p| [H3],[II],[5] [12], [IS]- They are all consistent with the requirement 
that the solutions of (1) and of SSE satisfy 

m = E[\m >< m\]- (2) 

where E[\ip(t) >< ip(t)\] is the expectation with respect to the distribution 
of the stochastic process \ip{t) >. The advantages of the description in terms 
of the pure states and SSE over the description by p are twofold. On the 
practical side, the computations are much more practical, as soon as the size 
of the Hilbert space is moderate or large [14]. On the theoretical side, the 
stochastic evolution of pure states provides valuable insides which can not 
be inferred from the density matrix approach [15] , [16] , [6] , [2] , [5] , [TTJ . 

There are two main approaches to the unravelling of the Linblad master 
equation: the method of quantum state diffusion [6J and the relative state 
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method [3J, 0, with specific advantages associated with each of the methods. 
The relative state method is usually used do describe the situations when the 
measurement is the dominant interaction with the environment. The method 
offers particular flexibility in that the master equation can be unravelled into 
different stochastic equations conditioned on the results of measurement. 
On the other hand the correspondence between the QSD equations and the 
Linblad master equations is unique, and is not related to a particular mea- 
surement scheme, or the form of the Markov environment. The resulting 
SSE is always of the form of a diffusion process on the Hilbert space of pure 
states, which is its main property to be explored in this paper. 

We shall concentrate on the unique unravelling of the master equation 
given by the quantum state diffusion equation, and explore the fact that it 
represents a diffusion process. QSD equation is the unique unravelling of (1) 
which preserves the norm of the state vector and has the same invariance as 
(1) under the unitary transformations of the environment operators {Li}, [6]. 
The equation is given by the following formula: 



\dip > = -iH\if> > dt 



+ 



£2 < i\ > U-i\U- < L\ >< L[ > 

I 



\^(t) > dt 



+ J2(Li-<L l >)\m>dW l (3) 
i 

where <> denotes the quantum expectation in the state \i/j(t) > and dWi are 
independent increments (indexed by /) of complex Wiener c-number processes 

The equation (3) represent a diffusion process on a complex vector space. 
We shall utilize the diffusion matrix of this process to define a Riemannian 
metric on the corresponding real space. We shall then study the properties of 
this diffusion metric as a field fixed by the environment and in relation to the 
stochastic evolution of the state vector, for different types of the environment. 
It will be shown, using examples of the dephasing and thermal environments 
and the measurement of an observable, that the curvature of the diffusion 
metric is a good indicator of the properties of the environment dominated 
evolution and its stability. We shall see that the curvature maxima of the 
diffusion metric coincide with the states that are preferred by the particular 
type of the environment. Furthermore, if the maxima are sharp and positive 
the stochastic dynamics governed by the environment and a Hamiltonian 
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perturbation that does not commute with Li, is likely to be attracted to the 
state with the maximal (positive) curvature. On the other hand the states 
that correspond to the negative values of the curvature are unstable. Our 
analyzes of the QSD equation, and the results, are strictly related to the 
fact that the equation represent a norm-preserving diffusion process, and in 
this sense are applicable to the stochastic modifications of the Schroedinger 
equation that describe a norm-preserving diffusion on the Hilbert space of 
pure states, like the QSD equation and, for example, the equations of the 
spontaneous collapse models [T3j . 

The structure of the paper is as follows. We shall first discuss, in the next 
section, a way to relate a Riemannian metric on a real space R 2n to a complex 
diffusion process on C n . Then, in section 3, we shall apply this procedure to 
define the Riemannian metric associated with QSD, and than the properties 
of this metric for various types of environments will be studied. Finally, in 
section 4, we shall summarize and discuss our results. 

2 Riemannian metric of a complex diffusion 

Using the following notation 

>) = -iH\1>> (4) 



+ 



IV >, 



£2 < ^|Lf|V >U- L\U- < ML\\^j >< V|^|V > 
. i 

B(\tP>)dW = Y,CLi-<^\Li\^>M>dW l . (5) 
i 



the QSD equation (3) assumes the standard form of a stochastic differen- 
tial equation (SDE) for an n-dimensional autonomous (stationary) complex 
diffusion process: 

d\i(j >= f(\iP >)dt + B(\i(j >)dW. (6) 

\ip{t) > and f(\ip(t) >) are complex vectors of complex dimension n, and dW 
are differential increments of an m-dimensional complex Wiener process: 

E[dW t ] = E[dWidW v ] = 0, 

E[dWidW v ] = 5 u ,dt, 

I = 1,2... m, (7) 
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where E[-] denotes the expectation with respect to the probability distribu- 
tion given by the (m-dimensional) process W, and W\ is the complex conju- 
gate of W\. B{\ip >) is n x m matrix, where m is at most n 2 — 1, and the 
diffusion matrix is 

G = BBK (8) 

Thus, G(\ip >) is Hermitian and nonnegative-definite. Notice that, unlike the 
case of a general SDE, the dissipative part of the drift (4) and the diffusion 
term (5) are determined by the same operators Li, and related in such a way 
that the diffusion equation preserves the norm of the state vector. 

The complex n-dimensional equation (3) generates 2n-dimensional real 
diffusion. Let us introduce the following real n dimensional vectors 

i 1 - 

p = 71^ 9=y|W + ^) 

t\) = -^=(q + ip), i) = -L(g-ip), (9) 

and a In dimensional vector X = (q,p). Similarly, we introduce real and 
imaginary parts of the vector / and order them as components of a 2n real 
vector T = (f R ,f I ), and introduce real and imaginary parts of the incre- 
ments of the complex m-dim Wiener process dW by 

dWi= (dW^ + idWD/V^, i = 1,2,... m (10) 

It is easily checked that the real and the imaginary parts are increments 
of a real 2m-dimensional process, i.e. 

EidW^dW*) = E(dW!dWj) = 5 iyj dt, E(dW t R dWj) = 0. (11) 

With this notation we have 



fdq\_ 1 / dip + # ' 
\dp) ~ v/2 Kdip-dipJ ' 



(12) 



Substitution of the complex equation (3) and its complex conjugate, leads to 
the following 2n dimensional real SDE: 

dq\_ff R (p,g)\ 1 (B R -B T \fdW R \ 
dp)~ \f I (p,q)) + V2KB 1 B R ){dW 1 )' [l6) 
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The matrix B of dimension 2n x 2m 

1 (B R -B J \ 

where 

(B) lJ = (B R ) ll +t(B I ) ij (15) 

gives the diffusion matrix Q for the real 2n dimensional diffusion described 
by the process (13) 

G-Btf- 1 ( (B R )(B*f + (Bi)(BY (B«)(BY-(Bi)(B«f 
^ ~ 2 V(S / )(S Jl ) T -(S fl )(S / ) T (5 7 )0BT + {B R ){B R ) T 

(16) 

We can write the matrix £ in terms of real and imaginary components of 
the n x n complex matrix G = BB^ as follows 



1 / G R G 1 \ 1 



g ~2{-G I G R J~2 



( G R \ , / -\ \ ( G 1 



V G R J \1 A G 1 



(17) 



where — G 1 = (G 7 ) T , since the matrix G is Hermitian. Furthermore, one 
can see that, besides the equalities between the entries corresponding to the 
symmetry of the matrix, there are other equalities 

(G)i,j = (G)i+n,j+n, i,j = 1,2. ..n (18) 

The matrix Q is symmetric and nonnegative, but it could be singular. 
However, the matrix Diag{l/2, 1/2, . . . , 1/2} + Q gives a Riemannian metric 
on the real 2n dimensional vector space. The factor 1/2 of the Euclidian 
part is chosen in order that the Euclidian norm of a vector corresponding to 
a complex n- vector of unit norm is also unity. 

Once the diffusion metric Diag{l/2, 1/2, . . . , 1/2} + Q is calculated the 
standard formulas [18] give the connection coefficients r* of the Levi-Civita 
connection for this metric in terms of the coefficients g^ = 5^/2 + (G)^ 
only 

= \g k \d,gx„ + d u g Xfl - d x g, u ) (19) 

Curvature tensor, Ricci tensor and the scalar curvature of the diffusion metric 
are also given by the standard formulas |18j : 

Rxi*v = d^vx ~ dv^nx + r^ A r^ — r^ A r^, (20) 
? A 



Ric^ = n = g^Bic^. (21) 
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Before we present the results of calculations of the diffusion metric and 
its curvature for different types of environments, we would like to consider 
briefly real representation of the QSD equation in the case when the Linblad 
operators are Hermitian. This includes, for example, the dephasing environ- 
ment or measurement, or the primary QSD [19] , [6] and other fundamental 
stochastic modifications of the Schroedinger equation [H] , [ID] , [12] ■ The goal 
of this digression is to point out to the connection between the general QSD 
equation (3) and some other stochastic modifications of the Schroedinger 
equation that have the form of a norm-preserving diffusion equation, and 
that consequently the construction of the diffusion metrics and its proper- 
ties are applicable to these equations also. In the case of Hermitian Linblad 
operators the real representation of (3) assumes a specially simple and illumi- 
nating form. Applying the same derivation as from equation (9) to equation 
(13) one obtains the following: 

dpi = —Hijqjdt + (2 < L > Lij — (L 2 )ij — < L 2 > 5ij)pjdt 

+ -j^ L H- < L > hj)PjdW R + ^(Lij- <L> 5 lJ )q j dW I , (22) 

where we have, for reasons of simplicity, included only one Linblad operator 
and the summation over repeated indexes is assumed. Noticing that for an 
arbitrary linear operator B 



BijQj — Oij 



d < B > 



BijPi = 5; 



d<B> 



i.i 



dpi 



(23) 



equation (23) becomes 



d < H > , dA 2 L , 1 
dpi = dt + ^ — dt + ^K 

oqt dpi V2 



d < L > rTTrR d < L > IT , 
dW R + — dW 1 



dqi 



dpi 



(24) 



where A 2 L =< L 2 > — < L > 2 . There is an analogous equation for dqi. The 
two sets of equations represent a diffusion process on R 2n , consisting of the 
drift given by a Hamiltonian dynamical system on R 2n with the Hamilton's 
function < H > and the dissipative part determined by A 2 L =< L 2 > — < 
L > 2 and the diffusion term determined by < L >. The drift and the diffusion 
are such that the norm of the vectors in R 2n is preserved. Furthermore, the 
equations are invariant under a global gauge transformation corresponding 
to the multiplication of vectors \ip > by a phase factor. Takeing into the 
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account the norm invariance and the global phase symmetry the equations 
can be written as a diffusion equation on the phase space S 2 ^ 1 / S 1 of the 
following form 

dX = QV < H > dt + V(A 2 L)dt + — |=V < L > dW (25) 

v2 

where V and QV are the gradient and the skew gradient on S 211 ' 1 / S 1 , and X 
denotes the set of 2n — 2 coordinates on the reduced phase space S 2n ~ 1 /S 1 . 
Equations like (25) have been analyzed as candidates for a description of the 
spontaneous state reduction in [13], or in the case L = H in [12J. 



3 QSD metric and qualitative properties of 
dynamics 

Application of formula (16) gives for the case (5) of the QSD equation an 
explicit procedure for calculation of the diffusion metric coefficients p^, in 
terms of the coefficients of the Linblad operators and the coefficients of the 
state tj) > in some bases \ip >= J2iCi\i >• The components of the diffusion 
matrix G = BB^ are given by 

m n 

B w (c,c) = EE( L " < L ' > l)fci(^ f - < ^ > Vk'fC^ (26) 

l 3,3' 

where: < L l >= J2 SS ' L ss 'c' s c s . Expressing q in terms of X\ . . . x 2n 

Xi = (q + q) / v2 i = 1, . . . n 
%i = \/— l(cj — Cj) / v^2 i — n + 1, . . . 2n, (27) 

separating of G R and G 1 and substituting in (16) finally gives the 4n 2 entries 
of the real matrix Q. 

We shall study the diffusion metric for the following three types of envi- 
ronments: (a) dephasing environment; (b) the environment corresponding to 
measurement of an observable and (c) thermal environment. The first two 
are represented by Hermitian and the third one by a non-Hermitian Linblad 
operators. The main geometrical object which we shall study are the diffusion 
metric norm of a state vector and its scalar curvature. In order to illustrate 
how these objects depend on the environment we shall use the simplest but 
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important quantum system, namely a single qubit. The system operators 
can be expressed as combinations of the Pauli sigma matrices a x ,a y ,a z , a 
state \ip > of unit norm is determined by < a x >, < a y >, < a z > or by the 
spherical angles (9, 4>) given by 

< & z > = cos(#) 

< a x > = sin(#) cos(0) 

< a y > = sin(#)sin(0), (28) 
The environment operators are |20j.[3] 

L = /ia + a- (29) 

for the dephasing and 

L = fi 1 a + + ii 2 <J- (30) 

for the thermal environment, with jMi and fi 2 proportional to the temperature, 
and finally for the measurement of, say, a z the Linblad operator is just 

L = ficr z . (31) 

The formulas for the entries gij of the diffusion metrics in terms of the 
coordinates X\, x 2 , x 3 , x A in the three considered cases can be conveniently 
written using the following notation: 

r2 2 I 2 j2 2,2 

— X I *^3? 2 — *^2 I X^ 

s = x\x 2 + X3X4, a = x\Xi — x 2 x%. (32) 

Because many of the metric entries are repeated, it is more convenient to 
present them in a list rather than to write down the corresponding matrices. 
Using the notation (32), the entries of the metrics in the three considered 
cases are: For dephasing: 

gu = l/2 + (/i 2 /16K(2 + ^)2 ! gi2 = (^2/16)^(2 + ^), 13 = O 

g 14 = (/i 2 /16)ad 2 (2 + d 2 ), 322 = 1/2 + 0* a /16)c#& g 23 = -g u 

924 = 0, 5-33 = g u , g M = g 21 , g AA = g 22 ; (33) 

For the thermal environment: 

gu = l/2 + d 2 [d 2 / u 2 + (2 + rf 2 ) 2 /i 2 ]/16, 

g 12 = s[dl(2 + d 2 2 )n 1 + (2 + dl)d 2 2 V2}/16,gi, 3 = 

9i4 = a[dj(2 + dj)^ + (2 + d 2 )4u 2 ]/16, 

g 22 = l/2 + rf 2 [d 2 /i 2 + (2 + rf 2 ) 2 /i 2 ]/16 

913 = 0, 5-23 = -9U, 924 = 0, #33 = 011, 934 = 921, 944 = 922, (34) 
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(a) (b) 




Figure 1: Poincare sections for the separability constrained non-symmetric 
quantum dynamics (39). The parameters are u> — 1, h — 1.5 and (a) /i = 1.3, 
(b) /i = 1.7 

and for the measurement of a z 

g n = l/2 + (g 2 /lQ)dl 2 (2 + d 2 1 -d 2 2 ) 2 , 

g 12 = (g 2 /16)s{d 2 .-d 2 1 -2)(d 2 1 -d 2 2 + 2), g 13 = 

gu = aidj + dj + dldl), g 22 = l/2 + (g 2 /16)d 2 2 (2 + d 2 2 -d 2 1 ) 2 , 

923 = ~gu, 924 = 0, g 3 3 = 011, 934 = 921, 944 = 922- (35) 

These formulas are used to compute the diffusion metric norm and the 
scalar curvature as functions of the state parameters 9 and 0. We shall first 
consider the dependence of the stated properties of the diffusion metric on 
the type of the environment and the coupling strengths /x, 1^1,1^2 and then 
analyze the relation between these properties and the stochastic dynamics of 
the state vectors. 
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Figure 2: Poincare sections for the separability constrained non-symmetric 
quantum dynamics (39). The parameters are u> — 1, h — 1.5 and (a) ji = 1.3, 
(b) » = 1.7 

In Figures 1 and 2 we illustrate the diffusion metric norm and curva- 
ture considered as functions on the sphere of states fixed by the type of 
environment and the value of the corresponding coupling /x, Hi,H2- Consider 
first Figure 1. The first row (fig. (a),(c),(e)) represent the diffusion metric 
norm and the second row (fig. (b),(d),(f)) the curvature for the three types 
of the environments and for some typical fixed values of the corresponding 
coupling strengths. The curvature is not constant, and can be positive or 
negative depending on the state vector and on the coupling strength. The 
maxima of the curvature can be sharp like in the cases of the dephasing and 
measurement of a z . On the other hand, in the thermal case the maxima is 
surrounded by a large neighborhood of states with almost maximal value of 
the curvature. Thus, the curvature has a sharp maxima only at the states 
which are clearly favored by the environment. If there are no such states 
the curvature maximum differs very little from the neighboring values. The 
curvature minima are at the states that are like repellers for the environment 
dominated dynamics. 

Dependence of the curvature maxima and the norm on the coupling 
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strength is illustrated in Figure 2 for the dephasing and the thermal en- 
vironments. The most important information from these Figures is that in 
the dephasing and measurement ( not shown) cases there are clearly sharp 
values of corresponding coupling strength where the curvature maxima goes 
from negative to positive values. Also, we see that the curvature minima are 
negative for all values of the coupling strength. 

We shall now study the relation between the sign of the curvature maxima 
and a stability of the stochastic dynamics of the state vector. The relation 
will not be analyzed in a mathematically rigorous way using an appropriate 
notion of the stochastic stability and considering the evolution of the met- 
ric as a stochastic process governed by the process \i/j(t) >. Instead, our 
strategy is to compute the curvature along different sample paths and see 
if the path remains near the state corresponding to the curvature maxima. 
We do such computations for the evolution governed by the environment and 
an additional fixed small hamiltonian, and we pay special attention to the 
case when the Linblad operators and the hamiltonian do not commute. The 
computations are repeated for the values of the coupling to the environment 
slightly above and below the critical value when the curvature maxima is 
zero. If the Hamiltonian perturbation is zero the sample paths that started 
near a maximum of the curvature remain near this maximum. For very 
small added Hamiltonian part and for a fixed value of the coupling to the 
environment, the sample paths of the system could wonder away from the 
maximum or could remain near it. In the former case we shall say that the 
stochastic dynamics is unstable and in the later case it is stable. The relevant 
computations are illustrated in Figures 3 and 4. 

In the case of the dephasing environment (or the measurement of a z ), 
when the maxima of the curvature are sharply picked, Figures 3 clearly illus- 
trate that positive curvature maxima correspond to the stability and negative 
to instability in the above mentioned sense. On the other hand, in the ther- 
mal case, the dynamics is always unstable even if there is no Hamiltonian 
perturbation. This is illustrated in Figure 4. We can conclude that the dif- 
fusion metric curvature provides us with a clear picture of the qualitative 
properties of the system's dynamics under strong influence of the environ- 
ment. 

It is well known that if the Linblad operators are Hermitian and commute 
with the hamiltonian, than the attractors of the stochastic QSD dynamics 
are the common eigenstates of the Hamiltonian and the Linblad operators 
|6j,[13j. The curvature maxima coincide with the eigenstates of the Linblad 
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Figure 3: Poincare sections for the separability constrained non-symmetric 
quantum dynamics (39). The parameters are u> — 1, h — 1.5 and (a) ji = 1.3, 
(b) /i = 1.7 




t 



Figure 4: Poincare sections for the separability constrained non-symmetric 
quantum dynamics (39). The parameters are u> — 1, h — 1.5 and (a) ji = 1.3, 
(b) /i = 1.7 
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operators, and consequently with the eigenstates of the hamiltonian. The 
probability of convergence to one of the attractors is, in this case, determined 
solely by the distance of the initial state from the attractor eigenstate, that 
is by the quantum mechanical transition probability, and does not depend 
on the parameters of the Hamiltonian and stochastic terms. The sign of the 
curvature maxima has no effect on this probability This is the reason why 
we expected that the relevance of the sign of the curvature maxima on the 
stochastic stability is manifested if the Linblad operator and the Hamiltonian 
perturbation do not commute. This expectation is qualitatively confirmed, 
as we described and illustrated in Figures 3 and 4, by numerical computa- 
tions. Observations of numerical sample paths, when the Linblad and the 
Hamiltonian operators do not commute, are enough to establish the quali- 
tative connection between the maxima of the curvature and the stability of 
small domains near the maxima. 

Finally, our treatment of the relation between the geometry of the dif- 
fusion and the stability of the stochastic dynamics is rather heuristic. We 
treated the diffusion metric as a given field on R 2n (determined by the Lin- 
blad operators), and we numerically studied the paths of the stochastic pro- 
cess \ip{t) > in relation to the sign of the curvature maxima. However, the 
problem of stability versus the properties of the diffusion metric should be 
formulated and studied using the appropriate notions of stochastic stability 
[2T],[22]. Nevertheless, we think that the numerical evidence strongly indi- 
cates that there is a clear relation between the sign and the shape of the 
curvature maxima and the systems dynamical stability. 

4 Summary and discussion 

According to the view of QSD theory, evolution of a state of an open quan- 
tum system is a diffusion process governed by a complex stochastic differen- 
tial equation on the Hilbert space of the system. The diffusion term of the 
QSD evolution equation explicitly depends on the operators modelling the 
environment and on the current state vector of the system. We have stud- 
ied the Riemannian metric associated with the diffusion term in the QSD 
equation. The metric is defined on the real 2n-dimensional space (here n is 
the complex dimension of the Hilbert space) and is directly related to the 
properties of the Linblad operators of the environment. We have shown that 
the scalar curvature of the metric has local maxima at states that are favored 
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by the corresponding environment. The curvature at different points, and 
in particular its local maxima, can be be negative or positive depending on 
the strength of the coupling to the environment. Also, the sharpness of the 
curvature maxima reflects the type of the environment. We have shown that 
there is a sense in which the sign of the curvature maxima is related to the 
stability of the corresponding state under the addition of a small perturba- 
tion that does not commute with the considered Linblad operator. If the 
environment type and the coupling strength are such that the curvature has 
sharp positive maxima, than the corresponding state is likely to attract the 
states of the system whose evolution is governed by the environment and a 
Hamiltonian that do not necessarily commute. On the other hand, if the 
curvature maxima are negative, the corresponding states are dynamically 
unstable under a small Hamiltonian perturbation that does not commute 
with the Linblad operators. In conclusion, the curvature of the diffusion 
metric is a relatively easy to calculate, and a very good indicator of what the 
environment dominated dynamics of the system would look like. 

The QSD equation describes the evolution of a pure quantum state using 
the Hilbert space of the quantum system, but, because it is norm-preserving, 
it gives also an equation on the state space, namely on the space of rays of 
the Hilbert space. The Riemannian metric associated with the diffusion on 
R 2n gives a Hermitian modification of the Fubini-Study metric on CP n ~ x . It 
is common to consider the complex projective manifold with the associated 
Fubini-Study metric as the proper framework for the geometry of quantum 
states [23], [21], so the modification of the metric due to the diffusion should 
also be formulated within this framework. 

The examples that we have analyzed in this paper are restricted on a 
single qubit under the influence of various types of environments. It would 
be interesting to analyze the properties of the diffusion metric in the case of 
coupled gubits, and in particular to see what is the curvature at the entangled 
states. Probably the proper framework for such analyzes is the formulation 
on CP™ -1 , mentioned in the previous paragraph, because the entangled states 
then have characteristic geometric interpretation [23J. 

Acknowledgements This work is partly supported by the Serbian Min- 
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FIGURE CAPTIONS 

Figure 1 Diffusion metric norm (a,c,e) and curvature (b,d,f) as functions 
of state parametrized by (9,(f>), for dephasing environment with fx = 0.6 
(a,b); measurement of a z with /j, — 1. (c,d) and thermal environment with 
[i x = 2,^2 = 1 (e,f). 

Figure 2 Diffusion metric norm (a,c) and curvature (b,d) as functions of 
9 for different values of the parameters fi or fii, fi2, and the maximum over 
(9, 4>) of the curvature as a function of // (e) or fii = 2/i2 = 2/x(/) . Figures 
a,b,e corespond to the dephasing and c, d, / to the thermal environment. 

Figure 3 Diffusion metric curvature (a,c) along the corresponding stochas- 
tic path illustrated in b,d for the dephasing environment and /i = 0.3 when 
maxi? < (a,b), and ji = 0.5 when max/? > 0(c,d). The small Hamiltonian 
perturbation is 0.01a x . 

Figure 4 Diffusion metric curvature (a) along a stochastic path (b) for 
the thermal environment and fii = 2/i2 = 1.6. The Hamiltonian part is zero. 
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